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A universal dimer is subject to secondary collisions with atoms when formed in a cloud of ultracold 
atoms via three-body recombination. We show that in a collisionally opaque medium, the value of 
the scattering length that results in the maximum number of secondary collisions may not correspond 
to the Efimov resonance at the atom-dimer threshold and thus can not be automatically associated 
with it. This result explains a number of controversies in recent experimental results on universal 
three-body states and supports the emerging evidence for the significant finite range corrections to 
the first excited Efimov energy level. 
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Few-body systems with a resonantly enhanced two- 
body scattering length a display universal properties in 
the sense that they are independent of the details of the 
short-range interaction potential [lj . In the two-body do- 
main with repulsive interactions, the most weakly bound 
energy level of the universal dimer depends solely on a, 
scaling as 1/a 2 . The central paradigm in the three-body 
domain, predicted in the early 70s by V. Efimov is 
associated with the infinite ladder of universal bound 
states with discrete scaling invariance. Beside a, univer- 
sal trimers depend on the three-body parameter, which 
accounts for all the short-range physics not already in- 
cluded in a. In recent experiments with ultracold atoms 
many aspects of universality were verified mainly by lo- 
calizing Efimov resonances which are associated with the 
crossing points of the trimers' binding energy levels with 
either the three-atom (for a < 0) or the atom-dimer (for 
a > 0) continua (for a recent review see Ref. Q). Surpris- 
ingly, positions of the lowest Efimov resonances for a < 
were found to be universally related to the van der Waals 
lengths of the two-body interaction potentials in differ- 
ent atomic species [4H14| indicating a universal three- 
body parameter. This remarkable experimental discov- 
ery has attracted intense theoretical attention which sug- 
gests a whole new understanding of Efimov physics in 
ultracold atoms 15, 0|. However, the situation with 
the Efimov resonances at the atom-dimer threshold (for 
a > 0) remains unclear. It appears to significantly devi- 
ate from the universal picture in all species investigated 
to date @, 0, [l3|, 17-20]. Suggestions have been made 
that the finite range of interaction potentials might be re- 
sponsible for these discrepancies. However, here we show 
that in some cases a different effect, directly related to 
the specifics of the experimental observables, can cause 
a significant shift in the resonances' positions. 

To show this effect we start by elaborating on two ex- 
perimental strategies developed recently to localize an 
Efimov resonance at the atom-dimer threshold, denoted 
in the following as a*. Both are based on the theoreti- 
cal prediction that resonant enhancements of the atom- 



dimer elastic (a e ) and inelastic (<jj) cross sections are 
expected in the vicinity of a* , according to the following 
analytical expressions [l], [2l[ : 



84.9 
20.3 



sin 2 [sq ln(a/a*) + 0.97] + sinh 2 77* 2 



sin' [sq ln(a/a*)] + sinh 2 77* 



sinh(2r7*) 



Ha 



3v sin 2 [s ln(a/o*)] + sinh 2 77* m A 



(2) 



where sq = 1.00624, 77* is the lifetime of the Efimov state, 
V is the velocity of the dimer and rriA is the mass of the 
atom [23]. In the first strategy, realized with ultracold 
Cs and 6 Li, an atom-dimer mixture is initially prepared 
in a shallow optical trap and the decay rate of dimcrs 
is then monitored as a function of a [la, [l^ - fl9l |. Thus, 
in this approach, the maximum of the enhancement in 
the atom-dimer inelastic collision rate is expected to co- 
incide with a*. In the second stra tegy , pioneered with 
39 K @ and then used with 7 Li 0, ill], the three-body 
recombination induced loss rate of atoms is measured. 
In each recombination event a universal dimer is formed 
with initial kinetic energy equal to one third of its bind- 
ing energy (Ed/ '3), which is usually much higher than 
the shallow optical trap depth. Moreover, the dimer is 
formed with the largest probability at the trap center 
where the density is highest. Then, on its way out, it 
may undergo secondary collisions with atoms, eject them 
out of the trap and cause the number of lost atoms per 
recombination event to greatly exceed three. Therefore, 
in this approach, the position of the maximum of the 
atom-dimer elastic collision rate is anticipated to reveal 
a ,„ . 

In this Letter we consider the problem of secondary 
collisions of dimers with atoms and develop a model to 
calculate their mean number as a function of a. Our 
central message is that the elastic and the inelastic pro- 
cesses are intimately interrelated and can not be consid- 
ered independently. According to Eqs. JT]) and ([2]), ai 
and <7 e are both enhanced when a = a*, however vari- 
ations in cr e may exceed those of Oi by orders of mag- 
nitude which impacts the experimental observables in a 
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nontrivial way. Appl ying this model to the experimental 
results of Refs. [y, |9|, l23j . we show that main controver- 
sies in the second experimental approach can be resolved 
when these processes are correctly included. We clar- 
ify the Efimov resonance' positions in 7 Li and 39 K and 
support the emerging evidence for significant finite range 
corrections to the first excited Efimov energy level in both 
systems 
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In the following we assume that the number of dimers 
is always very small compared to the number of trapped 
atoms and thus we neglect collisions between them. This 
situation is correct for all the experimental realizations 
considered here. We start with presenting two possible 
scenarios (Sc) that a dimer can undergo on its way out 
of the trap: 

Scl: The dimer collides elastically with k atoms and 
leaves the trap due to a high enough kinetic en- 
ergy. 

Sc2: After k clastic collisions, the dimer undergoes one 
inelastic collision, decays to a deeper bound state 
and leaves the trap without any further scattering 
events because o~ e and ai are enhanced only for the 
weakly bound universal dimers. 

To evaluate the probabilities of each Sc, we divide the 
total length / of the dimers' journey through the atomic 
cloud into N segments of infinitesimally small length 
81. Then, most generally, the probability to have an 
elastic/inelastic collision in the length segment 81 is: 
Pe/i = o'e/i^h where n is the mean density of atoms 28 1. 



Consequently, the probability to have no events in the 



same segment is: p n 



1 - Pi - p e 



According to Eq. t^fy, Oi depends on velocity and thus 
changes at each elastic collision because the dimer trans- 
fers a part of its initial kinetic energy to an atom. How- 
ever, it is instructive to start with the assumption that <7j 
remains constant independent of energy which makes the 
model clearer and analytically more tractable. Moreover, 
this simplification is reasonable in the limit of few colli- 
sional events and high initial energies of the dimer. Ex- 
perimentally it is realized in the region of relatively small 
scattering lengths and moderate atom densities which 
corresponds to some of the experimental conditions con- 
sidered later on. 

Energy independent (analytical) model. 

Scl: The probability for k elastic events to occur can 
be expressed by the choice of k segments out of total 
N, (,), in which the elastic scatterings occur with the 
probability p\ and multiplying it by the probability p^^, 
that in the remaining [N—k) segments no events happen. 
This probability is given by: 
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P {k; a e nl) exp (— OjnZ) . (3) 



where P (k; x) is the Poisson distribution. Note that 
although an inelastic event is not involved in Scl, the 
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FIG. 1: N(d) for the experimental parameters of Ref. (solid 
red line) and Ref. [t| (dashed-dotted line) . Parameters of the 
Efimov resonance are a, = 196a (marked with the dashed 
vertical line) and 77* = 0.1. 



probability pk,o does depend on m. This is because the 
probability of no inelastic event has to be included as a 
factor. 

Sc2: The inelastic event, which happens after k clastic 
collisions, has to occur at any of the N — k last segments 
with the probability pi. First, we fix the inelastic event 
at the m — th place where m £ [k + 1 , N] . Then, the 
procedure is similar to that of Scl. The probability for 
k elastic collisions to occur is expressed by the choice of 
k segments out of total m — 1, ( m Z )> with the probabil- 
ity p\ multiplied by the probability P™<^ n 1 e ~ k that in the 
remaining m — 1 — k segments no scattering events occur. 
Thus, the total probability for Sc2 is just the sum over 
all possible locations of m: 
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The mean number of scattering events in each Sc is 
then evaluated as the sum of weighted probabilities [29| : 
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fc=i 



N e ,i = 5^(fc+l)p fcl i = (l-exp(-0-iHl)) 

+ ^1 (1 - exp {-(Tiril) (1 + (70,)) . (6) 

Finally, the total mean number of events is then: N = 
N e ft + N e s- The central result of the paper is already 
seen in Eqs. (|5l6p : both of them include the interplay 
between elastic- and inelastic cross sections which strongly 
affects the experimental observablcs. 

Energy dependent model. Now we revamp the model 
taking into account the energy dependence of <Tj . When 
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a dimer is created via a three-body recombination, its 
energy is significantly larger than the kinetic energy of 
atoms. The mean energy of the dimer is reduced by a 
factor of 1/a 2 =5/9 per elastic collision with a nearly 



stationary atom |24( . As before, we calculate the proba- 



bilities and the mean number of events for each Sc. 

Scl: Let us consider k elastic events. Then the 
first one can occur only at any of the first N — (k — 1) 
segments with the probability p e . We fix its location 
at the segment j\ leaving the remaining ji — 1 seg- 



ments to have no scattering events with the probability 
(f — pice — p e ) Jl . Next, we fix the second scattering 
position j2 which then can only occur between segments 
ji + 1 and N — (k — 2) with the probability p e . In the 
remaining (j'2 — 1) — ji segments no events happen with 
the probability (l — pia 1 — p e )^ 2 n 1 and so on until 
the location jk is fixed. Finally, the total probability of k 
elastic events is evaluated by the sum of all these cases: 



N-k+l N-k+2 
Pkfi=p L e E E 



N k 

E IK 1 -^ 

Ji = l J2=ii+1 jk=jk-i+l m=l 



v Jm-Jm- 1-1 



(l - pia k - p e ) k , where j = 0. (7) 



Sc2: Here the probability is constructed almost iden- 
tically to that of Scl but with one variation. Before 
counting the elastic events we fix the segment where a 
single inelastic event occurs at the m — th position where 



m € [k + 1, N}. As a result, elastic events can happen 
only in the first m — 1 segments and we count them as in 
Scl. Thus, the total probability for Sc2 is just the sum 
over all possible locations of m: 
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where jo 
form: 



0. In the limit of 51 -> 0, cqs. (f7l8|) take the 



Pk.o 



k k 

Pk.i = — a 



3 



k ^ exp (- {criQp + <r e ) nl) 

h n™=o, m ^(« m -« j ) ' 

1 — exp (— (er e + aia? ) nl) 

o (fj+^')nLo, ro ^(« TO - ai ) 



(9) 
(10) 



which reduce to Eqs. (|3|4j) as a — > 1. The mean number 
of events N(a) can now be evaluated numerically based 
on the sum of weighted probabilities (see Eqs. (|5l6| )l. We 
now turn to apply the model to the recent experimental 
results which use the second experimental strategy. 

7 Li: Ref. [1| reports a* = 608ao, which should be cor- 
rected to a* ~ 410ao due to the shift in position of the 
Feshbach resonance [H, [22[. However, recent measure- 
ment of Machtey et.al. (MSGK) [23j], based on the same 
experimental technique reveals a* = 196(4)ao which is 
away from the Rice result by many times the experimen- 
tal uncertainty. There is, however, a major difference 
in the experimental conditions used by the two groups: 



while MSGK work with a low density thermal gas [H, [23| , 
the Rice experiment is performed on a Bose-Einstein con- 
densate (BEC). According to Ref. the mean density 
of atoms in the BEC is n « 5 x 10 12 cm~ 3 and half 
of the geometrical mean of the Thomas-Fermi radius is 
I p» 160 [im. In the MSGK case, for a typical temperature 
of T=1.4 yitK, radial and longitudinal trap frequencies of 
uj r = 2vr X 1.3 kHz and uj z = 2ir x 190 Hz and - 3.5 x 10 4 
atoms, the mean atom density is n w 1 x 10 12 cm -3 . The 
geometrical mean size of the atomic cloud is I s=s 9.5 /xm 
(one standard deviation of the Gaussian distribution). 
Therefore the column density (nl) in the Rice experiment 
is larger by a factor of ~ 80. 

For the second strategy, all Sc's contribute to the ex- 
perimental observable and so we show N(a) in Fig. Q] 
for both experiments. We used the energy dependent 
model with a* = 196a . The solid red (dashed-dotted 
blue) line corresponds to the MSGK (Rice) experimen- 
tal conditions. In the MSGK case, the maximum in N 
nearly perfectly coincides with the position of a*, which 
indicates that the experimental observable indeed reveals 
the position of the Efimov resonance. However, for the 
Rice experimental parameters, the maximum of TV is not 
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FIG. 2: N(a) for the experimental parameters of Ref. [6|. In- 
set shows a peak at low scattering length which also reason- 
ably agrees with the experiment. Parameters of the Efimov 
resonance are a, = 515ao ( marked with the dashed vertical 
line) and 77* = 0.1. 



at a* but rather appears at ~ 375ao, in very good agree- 
ment with the reported result Q. A simple intuitive 
understanding of this difference can be obtained when 
the limit of aifil « 1 (» 1), roughly matching MSGK 
(Rice) experimental conditions, is taken in Eqs. (|5I6[) . 
Then in the MSGK case, TV rj a e nl, which corresponds 
to the mean number of the elastic collisions and peaks 
at the maximum of er e , i.e., at a = a*. However, for the 
Rice case, N « 1 + <r e /ai, which is maximized at a dif- 
ferent position, namely when srjln(a/a*) + 0.97 = n/2 
(see Eqs. ([Tp]) ) (30l |. This corresponds to a/a* « 1.8 
in agreement with the energy dependent model and the 
Rice result. 

Note that the model predicts a larger width of the N(a) 
maximum for the Rice experimental conditions, qualita- 
tively well identified in the experimental results 0, [23| 
and, thus, further strengthening our model. However, 
quantitative comparison requires inclusion of additional 
conditions such as reduction of the dimcr energy below 
the trap depth after a certain number of clastic collisions. 
Here we only associate the maximum in the three-body 
recombination loss rate with the maximum in the mean 
number of secondary collisions. In the region of large a 
not all of them, if any, lead to atom loss, which cuts down 
the long tail of N(a). This is, however, beyond the scope 
of the present discussion. 

39 K : In Ref. 0] two Efimov resonances at the atom- 
dimer threshold are reported at a*i = 30.4ao and a*2 = 
930ao. We shall concentrate our analysis on the second 
resonance as the first one is deeply inside the nonuniver- 
sal region where a*i < ro with ro = 64.5ao being the 
van der Waals length of 39 K. The experiment is also per- 
formed in a BEC with n rj 6.5 x 10 12 cm~ 3 and half 
of the geometrical mean of the Thomas-Fermi radius is 
I w 7 /im. As before our model predicts a significant shift 
in the real position of the resonance as compared to the 



TABLE I: Comparison between newly extracted values of the 
Efimov resonances at the atom-dimer threshold (a*) and the 
values predicted by universal theory (a^ T ). While a^ 1 de- 
notes the originally reported experimental results, a, shows 
the corrected values of a* 1 according to our model if appli- 
cable, a*) indicates positions of recombination minima in the 
measured three-body recombination spectra and is used to 
find a^ T via a universal relation a^ T ~ l.la^ (22.7) n_1 ^ 2 
where n can be either or 1 for the purpose of this table [lj|. 
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- 410 
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515 
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la, 



397 = af 



1.48 
~ 1.4 
2.28 
2.53 
2.77 



32.5 

64.5 
101 



"See also ref. [S|]. 

6 The values are rescaled according to correct Feshbach resonance 
parameters, i.e. position and energy width [iil |22||. 

c See also refs. [j, 0. Note that is reasonably well related to an 
Efimov resonance at a < Q| despite the fact that it is too close 
to ro and thus might be subjected to the finite range corrections. 



reported one. To agree with the experimental value of 
~ 930ao we need to take a» = 515ao, which reflects the 
same numerical factor of 1.8 as before. N(a) is shown in 
Fig. [2] In the inset we represent the second peak located 
at ~ 23ao and calculated with the slightly different ex- 
perimental parameters in accordance to Ref. Q but with 
the same a*. Although it is close to the measured value it 
might be merely a coincidence. Note also the difference 
in widths of the two features which again qualitatively 
agrees with the experiment. 

We can now compare the newly extracted positions of 
the Efimov resonances at the atom-dimer threshold for 
both species with the values predicted by universal the- 
ory (see Table H| . In both cases a* is notably shifted to 
lower values. MSGK associate this shift with the mani- 
festation of finite range corrections 0, [25[ . For 7 Li the 
shift is smaller than for 39 K. In both cases the Fesh- 
bach resonances are of intermediate character between 
being closed and open channel dominated, but r of 7 Li 
is about half that of 39 K . It is interesting to note that 
in Cs the shift is even larger despite the open channel 
dominated character of the Feshbach resonance (see Ta- 
ble U) 0, UtI [Lsj . However, the ro of Cs is the largest 
among the other species and we can then conclude that 
the recent experimental results support a scaling of the 
finite range corrections in ro. Finally, we note that exper- 
iments provide growing evidence for a good agreement of 
the position of the lowest minimum in the three-body re- 
combination spectrum with universal theory despite the 
fact that it is measured at even lower scattering lengths 
than a* 1 which still remains a puzzling question. 
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